Order-by-disorder in the XY pyrochlore antiferromagnet revisited 



O 

bj:): 

< 

o 

cn 

o 

S 



I 

o 



m 
o 

00 
O 



X 



Pawel Stasiak,'-^ Paul A. McClarty,'-^ and Michel J. P. Gingras'-'* 

' Department of Physics and Astronomy, University of Waterloo, Waterloo, ON, N2L 3G1, Canada. 
^Department of Mathematics, University of Reading, Whiteknights, PO Box 220, Reading RG6 6AX, UK. 
^Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str 13, Dresden, 01187, Germany. 
''Canadian Institute for Advanced Research, 180 Dundas Street West, Suite 1400, Toronto, ON, MSG 1Z8, Canada. 

(Dated: September 1, 2011) 

We investigate the properties of the XY pyrochlore antiferromagnet with local (111) planar anisotropy. We 
find the ground states and show that the configurational ground state entropy is subextensive. By computing 
the free energy due to harmonic fluctuations and by carrying out Monte Carlo simulations, we confirm earlier 
work indicating that the model exhibits thermal order-by-disorder leading to low temperature long-range order 
consisting of discrete magnetic domains. We compute the spin wave spectrum and show that thermal and 
quantum fluctuations select the same magnetic structure. Using Monte Carlo simulations, we find that the state 
selected by thermal fluctuations in this XY pyrochlore antiferromagnet can survive the addition of sufficiently 
weak nearest-neighbor pseudo-dipolar interactions to the spin Hamiltonian. We discuss our results in relation to 
the Er2Ti207 pyrochlore antiferromagnet. 

PACS numbers: 75 . 1 0.Dg 75 . 1 0. Jm 75 .40.Cx 75 .40.Gb 



I. INTRODUCTION 

The geometric frustration of magnetic interactions on lat- 
tices of magnetic moments often leads to a configurational 
classical ground state entropy that scales with the volume of 
the system, V, as V"' with < a < 1. This can have some 
unusual consequences. A well known example is the Ising 
model on the triangular lattice with nearest-neighbor antifer- 
romagnetic interactions which has an extensive ground state 
entropy and exhibits no finite temperature transition. ' While 
in real materials a V" entropy left by the leading interactions 
is often energetically lifted by weaker interactions, leading to 
long-range magnetic order, there are some exceptions. For 
example, in the Dy2Ti207 and Ho2Ti207spin ice materials,- 
in which the magnetic moments are described by Ising spins, 
the extensive (a = 1) low temperature entropy caused by frus- 
tration of the leading effective ferromagnetic nearest-neighbor 
interactions is indeed lifted by the perturbing long-ranged part 
of the dipolar interaction!^ However, the degeneracy lifting 
in this system is so weak that the theoretically expected phase 
transition to magnetic long-range order is inhibited by a freez- 
ing into a spin ice state without long-range order^^ 

Another possibility in a system with an exponentially 
(cxp[Cy"]) large number of classical degenerate ground 
states is that thermal or quantum fluctuations might select a 
subset of states about which the density of zero modes is 
greatest. These entropic and quantum state selection mech- 
anisms are both referred to as order-by-disorderj^^ii Among 
pyrochlore antiferromagnets, in which the spins sit on a lat- 
tice of corner-sharing tetrahedra, Moessner and Chalker have 
given a criterion for the occurrence of long-range order in- 
duced by thermal fluctuations.'^ This criterion is based on 
the degree of divergence of the statistical weight of particu- 
lar spin configurations — a power-counting argument depend- 
ing on the number of zero energy excitations (zero modes) 
for a given spin configuration and the number of dimensions 
of the ground state manifold. For example, this criterion in- 
dicates that the XY antiferromagnet with globally coplanar 



spins (spins perpendicular to the global [001] axis) should ex- 
hibit entropic selection — a result which is borne out by Monte 
Carlo simulations!^ For such XY systems, this comes about 
because the number of zero modes about collinear spin con- 
figurations is proportional to the number of spins whereas the 
configurational entropy in the ground state is subextensive, 
growing as y2/3^nj3 

This article is concerned with the pyrochlore XY antiferro- 
magnet with local (111) spins meaning that there is a differ- 
ent easy plane for each of the four tetrahedral sublattices»i^"— 
Because such a model preserves the cubic symmetry of the 
pyrochlore lattice and because the single ion crystal field can, 
and does in various materialsjii generate such an anisotropy, 
it is more physical than the aforementioned pyrochlore XY 
model with a global easy axis.'^ The model has been rec- 
ognized to exhibit a continuous degeneracy in its classical 
ground stateJ^ii^ Monte Carlo simulations of the local (111) 
XY antiferromagne t ' ^- ' ^- ' ^ indicate that it exhibits two phases 
— a high temperature paramagnetic phase and a low tempera- 
ture long-range ordered phase. We refer to the magnetic struc- 
ture in the ordered phase as ^2 to be consistent with Ref. Iisll 
and the group theory literature. A calculation of the spec- 
trum of the Hessian about different discrete ground statesi^ 
suggests that the observed long-range ordered spin configura- 
tion in Monte Carlo simulations has the largest density of zero 
modes of all the degenerate ground states and, consequently, 
that the observed transition is an example of classical (ther- 
mal) OBDji^ii^ However, it has been suggested that selection 
of the long-range ordered ■02 state might not survive in the 
thermodynamic limitji^ 

In this article, we give a systematic account of the prop- 
erties of the local (111) XY pyrochlore. We present in Sec- 
tionlUthe model and its ground states. In Section HIH we dis- 
cuss some of the details of the Monte Carlo simulations per- 
formed in this work. Section|IV]reports results of an analytical 
and numerical investigation of the thermal order-by-disorder 
mechanism, providing strong evidence that the fluctuation se- 
lection mechanism of the 02 state does survive in the ther- 
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modynamic limit and giving further insight into its physical 
origin. We also include in Section [TVl a subsection showing 
that there is a quantum order-by-disorder mechanism in the 
XY model as speculated but not shown in Ref. lfl4|] . Finally, 
in Section |V] we describe the material Er2Ti207 which is an 
easy plane antiferromagnet exhibiting the ijj2 structure in its 
ordered phas o'"^i'^ and which provides an experimental moti- 
vation for studying this model. In particular, we discuss the 
effect of weak dipolar interaction on the XY antiferromagnet 
and the problem this interaction poses for understanding the 
long-range ordered phase of Er2Ti207 with ip2 structureJ^ii^ 



tetrahedron). Then, we label these branches in the following 
way: 



Branch 1 
Branch 2 

Branch 3 
Branch 4 



il = 
ii = 



-03 = --^4 



- 47r 



— V2 



(2) 



II. MODEL 

In this work, we mostly focus on the problem of the zero 
and finite temperature behavior of interacting classical spins 
of length |S| = 1 on the sites of a pyrochlore lattice of corner- 
sharing tetrahedra with an infinite single-ion anisotropy such 
that the spins lie within their respective local XY planes per- 
pendicular to the local (111) directions. In Section [TV B| we 
discuss the problem of order-by-disorder due to quantum fluc- 
tuations in a model with spin operators S 

The interactions are taken to be antiferromagnetic isotropic 
exchange between nearest neighbors with coupling J (J > 0). 
Later on we also consider, as a perturbation, pseudo-dipolar 
interactions solely between nearest neighbours and with cou- 
pling strength V. Thus the Hamiltonian is taken to be 



\ ^ Sj • Sj 3(Sj; • R.ij)(Sj • Ri 



(1) 

where R^^-^ is the nearest-neighbor distance. 

Consider first the exchange-only model with V — Q. In 
this case, the Hamiltonian in Eq. ([T]i can be put into the form 
Hex = JJ2t (^t ~ 45*^), where the sum runs over all con- 
nected tetrahedrai^ and St is the total spin on each tetrahe- 
dron. It follows that the ground states are all those states 
with zero net magnetic moment (St — 0) on each tetrahe- 
dron. Therefore, we write down the conditions for the three 
components of the total moment on a tetrahedron to be zero. 
In doing so, we impose the XY constraint so that the orien- 
tation of spin a, (for sublattices a = 1, 2, 3, 4), is given by a 
single angle (j)a measured with respect to axes within the local 
plane (normal to the relevant local [111] direction) given in 
Ref. 111911 . The condition of zero moment on each tetrahedron 
can then be written as 



cos (0i) + cos (02) ^ COS (03) + cos (^4) 
cos + COS (03) = cos (02) + cos (04) 

cos (0") + cos (0") = cos (02) + cos (03^ 



where 



There are four so- 



= 0, + f and 0:,' = 
lution branches to these equations. Each branch corresponds 
to a continuous degeneracy wherein all four spins are rotated 
smoothly within their respective local [111] XY plane. 

We place an overbar on 0a {(pa) to signify the angle for 
sublattice a giving an energy minimum (zero moment on each 



A further discussion of these solutions can be found in Ap- 
pendix A. 

To enumerate all the ground states on the pyrochlore lat- 
tice we first tile all the tetrahedra with a particular spin con- 
figuration from Branch 1. Then, we choose a line of nearest- 
neighbor spins traversing the length L of the system. The sub- 
lattice labels of the spins on the chain alternate between two 
values a and b. There are six such pairs of labels. One can then 
transform the spins along the chain so that the spin configu- 
rations of the associated tetrahedra belong to another branch 
of solutions. For example, consider a single chain made of 
sublattices 3 and 4. All the local angles along this chain are 
identical initially and equal to, say 6. We can transform these 
to —9 with no energy cost. Therefore, the entropy within the 
ground state manifold scales as as first noted in Ref. ifisll . 
This is in contrast to both the Heisenberg pyrochlore antifer- 
romagnet and the global easy axis (Ising) pyrochlore antifer- 
romagnet both of which have an extensive entropy. 

We note that the four branches in Eq. (|2]i intersect in pairs. 
These intersection points are at special sublattice angles = 
7i7r/3 with integer n. We refer to these as ?/'2 states in the rest 
of this article. By exploiting these intersection angles to move 
between the branches, one can smoothly visit all the ground 
states on a single tetrahedron and, indeed, on the whole py- 
rochlore lattice. If we return to the above chain of sublattices 
^3 and #4, the = configuration allows the tetrahedra 
along this chain to pass smoothly from Branch 1 to Branch 
2. As shown in Refs. Ulisll and in Section HV] below, ther- 
mal fluctuations have the effect of selecting a magnetic struc- 
ture with q = ordering wavevector and spin orientations at 
these discrete angles. There are six distinct ijj2 ground states 
which are the six q = ordered states with tetrahedra tiled 
with local angles 0a ~ mr/S for sublattices a = 1,2,3,4 
and with integer n. One can take the observation that the lat- 
tice zero modes are along sublattice chains to understand an 
aspect of the Monte Carlo results of Ref. ifisll ; in particular, 
the finite-size scaling of the average energy of the -02 states at 
low temperature. Since this point is somewhat removed from 
the main story of the paper, we present the argument in Ap- 
pendix [C] 

When 2? 7^ 0, the continuous ground state degeneracy of 
the exchange only model is "immediately" replaced with a 
discrete global degeneracy with q = ordering wavevector 
selected from the manifold of states described above. These 
energetically selected states are referred to as the iJj a sta tes — 
or Palmer-Chalker states in the literature after Ref. ll20ll . The 
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angles specifying the 1/^4 states are 

State 1 : = 02 = ^ 4>3 = 4>4 = 

- - Ttt - - it 

State 2 : (pi = (p4 = — 92 = 93 = -p 

6 6 

2^ Ivr Stt 

State 3 : = (jj^ = —— 02 = (/)4 = — (3) 

6 6 

and the time-reversed configurations. In anticipation of what 
follows in Section |Vl we note that the '04 states are the 
ground states one finds for antiferromagnetic nearest-neighbor 
exchange with sufficiently weak nearest-neighbor pseudo- 
dipolar as well as for true \/r^ long-range magnetostatic 
dipolar interactions in the classical Heisenberg pyrochlore 
antiferromagnet modeli^— Interestingly, the ^■'4 states are 
found experimentally to be the ground state of the Gd2Sn207 
pyrochlore antiferromagnet^! but not of the closely related 
Gd2Ti207 materiali^ 



III. MONTE CARLO METHOD 

In Sections flV Cl andFVlbelow. we report results from Monte 
Carlo simulations of the local (111) XY pyrochlore antiferro- 
magnet. In this section, we give details of the Monte Carlo al- 
gorithm and the observables that were measured in the Monte 
Carlo simulations. 

The Monte Carlo simulations were performed using par- 
allel tempering^- in which Nt replicas of a system of N 
spins, each at a different temperature and with a different 
series of pseudo-random numbers, are simulated simultane- 
ously. In addition to local spin moves, parallel tempering 
swaps that exchange configurations between a pair of tem- 
peratures are attempted. The configuration swap attempts are 
accepted or rejected based on a Metropolis condition that pre- 
serves detailed balance. Parallel tempering has been shown, 
in systems known to equilibrate slowly using other meth- 
ods, to improve performance substantially;^ Replica swaps 
are attempted with a frequency of one attempt every 100 local 
Monte Carlo sweeps. A local Monte Carlo sweep consists of 
N spin move attempts. In our simulations, Nt ~ 64 with 
either a constant increment between the temperatures, or with 
the temperatures self consistently adjusted to obtain a uniform 
parallel tempering acceptance rate. 

With each spin carrying a single angular coordinate 0^, the 
local spin moves involve choosing an angle increment 5(j)i 
from a uniform distribution between — 50niax and J0max- The 
angle of spin Sj was updated to 0j + (50; and each tentative 
spin rotation was accepted or rejected based on a Metropo- 
lis test. The maximum increment (50max was updated every 
100 Monte Carlo moves in order to maintain the spin move 
acceptance rate at 50%. 

Physical observables were computed every 100 Monte 
Carlo sweeps. To determine the presence of long-range or- 
der with ordering wavevector q = (expected for sufficiently 
small 2?/ J in the model discussed above2£i22), the sublattice 
magnetization was computed^! 




where each spin carries an fee lattice label i and a sublattice 
label a (see Ref. H^l) and the number of sites in the lattice is 
N = ANp, where Np is the number of fee sites. The angled 
brackets (. . .)th denote a thermal average. In order to distin- 
guish the i/)2 (Refs. ITilsIl ) and ^4 (Refs. HMl) phases, 
we introduce unit vectors e'a'^'^^'^ which are oriented in the ex- 
pected spin directions on each sublattice a for magnetic struc- 
ture identified by the label 7 with the domains labelled rf, for 
both the 7 = "02 and 7 = -04 structures. From the combina- 
tion 

Np 4 

^(.(.^^.^^^S^.-ei-C^)), (5) 

P i=l a=l 

we compute the order parameter 




for the = tp2 (Refs. HTs!]) and 7 = -04 (Refs. fT8'.'2^) 
magnetic structures. In Eq. (|6]l, the sum is taken over a choice 
of three out of the six magnetic domains, for each of these 
two structures, which are not related to one another by time 
reversal. The spin directions corresponding to the domains 
for 014 and V'2 are given in Section II. — for 04 in Eq. Q and 
for 0)2 we have all angles 0^ = mr/3. The order parameter 
for -0 4, is the same one computed in the simulations of 
Ref. 1I27I1 . The magnetic specific heat per spin was computed 
from the fluctuations in the total energy of the system. 

The sensitivity of the results to the initial spin configura- 
tions was assessed by comparing the results of simulations 
starting from (i) random configurations with a different con- 
figuration for each thermal replica, (ii) 0^2 ordered states 
and (iii) V'4 ordered states. To ensure that equilibration was 
reached for each simulation, we checked that the results were 
independent of initial conditions. Also, the evolution of the or- 
der parameters was monitored during the course of each simu- 
lation to ensure that they reached a stationary state before the 
statistics were collected. Equilibration issues are discussed 
further in Sections HVCl andlVl 

IV. ORDER-BY-DISORDER 

In this section, we consider the exchange-only model {T> = 
0) given in Eq. ([T]i. General arguments given in Ref. lfl2ll indi- 
cate that the XY antiferromagnet with coplanar spins should 
exhibit a thermally driven order-by-disorder transition. This 
argument does not straightforwardly carry over to the non- 
coplanar (111) XY antiferromagnet. Simulation evidence for 
thermal order-by-disorder in the local (111) XY model transi- 
tion was presented in Refs. lfl3l - [l6ll . The classical degenera- 
cies of this model were identified in Ref. ifTsIl and order-by- 
disorder was found via Monte Carlo simulations. However, 
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the possibility was mentioned in Ref. fl55| that the tempera- 
ture at which long-range order with a nonzero ^'2 order pa- 
rameter develops might vanish in the thermodynamic limit. 
We present simulation results which provide compelling evi- 
dence that, for the exchange-only (2? = 0) model, a first order 
phase transition to a long-range ordered ip2 state persists in 
the thermodynamic limit. We begin, however, with a previ- 
ously unreported calculation of the free energy including only 
harmonic fluctuations which exposes a thermal ?/'2 order-by- 
disorder in the thermodynamic limit. Then, having investi- 
gated the order-by-disorder mechanism in the exchange-only 
{T> = 0) model, we discuss in SectionlVlthe effect of compet- 
ing nearest-neighbor pseudo-dipolar (T) 7^ 0) interactions in 
this model. 



A. Computation of the free energy 

In this section, we show that certain discrete spin configu- 
rations from the manifold of q = ground states minimize 
the free energy computed from harmonic fluctuations about 
the classical ground states. We assume that every tetrahedron 
on the lattice is tiled with the same spin configuration (i.e. 
that the ordering wavevector is q = 0). If we constrain the 
ordering to be q = 0, the spin configuration is fixed by spec- 
ifying four angles (f>a - one for each sublattice. Let the an- 
gles in a ground state configuration be denoted 0a for which 
the ground state energy is H{4>a) = NEg, where N is the 
number of spins. We then we consider small fluctuations Scpi 
about these angles 0^ = (f)i + (50^. The terms linear in 5 (pi van- 



ish, so the Hamiltonian H = NEq + H2 



where H2 is 



the part harmonic in the angular deformations. H2 is written 
in k space as H2 = Ek,a,6 '5'^a(k)A«^(k)50b(-k) . Here, 

<50a(k) = (1/ViVp) E'j^'^ CXp(ik • (R^ + ra))50a(Rp), 

where are the fee lattice points and Ta are the vectors for 



the tetrahedral basis (see Ref. 112611 for notation convention). 
This choice of convention for the lattice labelling ensures that 
the Hessian A"''(k) is real. The eigenvalues AA(k) of A''''(k) 
are nonnegative, reflecting the stability of the ground states. 
The spectrum of A°'^ is computed in k space as a function of 
the ground state for each branch. 

One finds that for the special minimum energy configura- 
tions 0a = n7r/3 for a = 1,2,3,4, the four eigenvalues, 
AA(k), {A = 1, 2, 3, 4) of the A"''(k) Hessian take the form 

0a = 0, TT Aa = 1 ± cos(k • ri2), 1 ± cos(k • r34) 

0a = 7r/3, 47r/3 A^i = 1 ± cos(k • ria), 1 ± cos(k • r24) 
0a = 27r/3, 57r/3 A^ = 1 ± cos(k • r23), 1 ± cos(k • ri4) 

where each row gives the four eigenvalues for the indicated 
particular set of 0a (a — 1, 2, 3, 4) angles and which corre- 
spond to the aforementioned 11^2 states. The vector Yah joins 
nearest neighbors with sublattice labels a and 6. The ■4>2 states 
are distinguished from the other ground states in having a 
much higher density of zero modes — 2 planes of zero modes 
in the first Brillouin zone.'^-'^ The planes come about because, 
at these angles, one can smoothly introduce defects into the 
system along (what are usually referred to as a and /3) chains 
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no. 1 : (Color online) Plot showing the harmonic free energy con- 
tribution I2 for each of the four branches {A = 1, 2, 3, 4) of ground 
states with the 0a angle for each sublattice a = 2, 3, 4 expressed 
in terms of 0i via the parametrization given in Eq. Each of the 
four curves is labelled by its ground state branch number taken from 
Eq. (|2}. Note therefore that the horizontal axis does not identify a 
unique spin configuration but rather a unique configuration for each 
of the four branches. The minima in the free energy appear where 
pairs of ground state branches meet - at the ^2 spin configurations. 



in the crystal as described in Section|II] At angles away from 
the V^2 states, one can still introduce chain defects but not con- 
tinuously. 

We now compute the free energy for each q = configu- 
ration after dropping all terms in the Hamiltonian beyond the 
harmonic terms. The free energy F[0] at inverse temperature, 
j3, is given at the harmonic level by 



F[0] 



NEg - i log 



n j d[<50a(k)] j exp(-/3iJ2) 



(7) 



and hence 



F[0] « NEg 



N fn 



^^log(dctA(k)). (8) 



In the limit as N 
right hand side of 



Eq. 



00, 

m 



the 

is 



third 

N_ 1 

4 



(last) 



term on 
with I2 



the 



4 (27I-)3' 

J2aI '^'^klog (AA(k)) for the four ground state branches 
given in Section |ll] (see Appendix iBl for a discussion of the 
various prefactors of 12). We evaluate the integral numeri- 
cally using a Monte Carlo method with 10* points, noting that 
the singularities for the ip2 spin configurations are integrable 
because the integrals take the form dk log k for some con- 
stant c.^ The results are shown in Fig.[T] Evidently, I2, and 
consequently the free energy F[4>], is minimized for the i/;2 
states at 0a — mr/S where pairs of branches meet. 

It is often the case that one can simulate the effect of order- 
by-disorder by introducing a term into the Hamiltonian of the 
form 



OBD 



IE 

2, J 



Si ■ So 



(9) 



5 



that selects the most collinear spin configurations among the 
classical ground states (see, for example, Refs. lIsUloll). The 
usual argument for the selection of such states is that colhnear 
spin configurations have, among all states, fluctuations that 
couple most strongly because fluctuations are responsible for 
effective fields perpendicular to the spin direction even in the 
broken symmetry phase. The local XY and zero moment con- 
straints of the XY pyrochlore antiferromagnet ensure that the 
spins cannot be collinear, but it is interesting to ask whether 
the ijj2 configurations are the most collinear states within the 
set of ground states. One finds that Eq. (|9|l is constant within 
the whole ground state manifold of Eq. (|2]i. However, the sum 
— J2i j l^i ' Sjl is minimized by the ?/'2 states which lends 
some credence to the intuition that the most collinear states 
among all the classically degenerate zero temperature ground 
states must be selected. 
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B. Quantum selection 

Having shown that thermal fluctuations select the ijj2 states, 
we now turn to the effect of quantum fluctuations which, in 
general, need not select the same states. In this section, we 
present the spin wave spectrum computed using the Holstein- 
Primakoff transformation treated in a large S expansion and 
truncated at harmonic order. The calculation is performed 
for the Hamiltonian in Eq. ([T]i with V ~ 0. In the coordi- 
nate system with z axes taken along the (111) directions, the 
Hamiltonian is written H = ^J^^SfS^ where a and /3 
denote the spin components. The local Ising components of 
the matrix of interactions J"-^ are set equal to zero — this 
imposes a soft XY constraint because the computation of the 
spin wave spectrum implicitly allows fluctuations out of the 
easy planes. This is relevant to the Er2Ti207 pyrochlore an- 
tiferromagnet whose single-ion crystal field doublet, charac- 
terized by an anisotropic g-tensor with two eigenvalues such 
that g±^ > , allows for a description in terms of an effective 
spin- 1/2 modelj22, Working in reciprocal space with N spins 
and Np primitive lattice sites, one rewrites the spin Hamilto- 
nian in terms of boson operators, with2L2^ 




on each site with the new z axis now taken to be the quanti- 
zation axis. The quantization axis is taken within the ground 
state manifold of the model with antiferromagnetic exchange 
as parametrized by local angles (f>a for a ~ 1,2,3,4 given in 
Eq. (|2]i. There are four flavors of bosons corresponding to the 
distinct sublattices labelled with subscript a. One performs 
a Bogoliubov transformation taking boson operators ajj(k) 
and tta (k) into spin wave creation and annihilation operators 



FIG. 2: (Color online) Zero point energy for the pyrochlore XY 
model from the linear spin wave spectrum computed from the clas- 
sical ground states. The horizontal axis is the (j> parameter given 
in Eq. ^ and each curve carries a label identifying the branch of 
ground states to which it belongs. The quantum correction to the 
classical ground state energy is minimized for the tp2 states. 

(k) and ca (k) so that the Hamiltonian to harmonic order is 
brought to the form 

H[^] ^ ~NJS{S + 1) + JS^eAO<) 

+ JS'^eA(k)c^(k)cA(k) (13) 

where eA(k) are the spin wave energies. Further details of 
Holstein-Primakoff linear spin waves on a pyrochlore lattice 
of spins can be found in Ref . ll2Tll . As one would expect for 
an antiferromagnet, the dispersion for the model Eq. ([T]i about 
the zero modes is linear in |k|. Just as in the classical case, 
the zero modes appear in pairs of planes in the first Brillouin 
zone for the ip2 states. The harmonic correction to the ground 
state energy is 

which we have evaluated numerically. The results are shown 
in Fig. |2] To harmonic order, one observes that, among the 
q = ground states, the zero point energy is minimized at 
the ip2 spin configurations {(j>a = titt/S for all a), so a quan- 
tum order-by-disorder mechanism selects the same states as 
thermal fluctuations. 



C. Monte Carlo results 

To confirm the thermal order-by-disorder mechanism ar- 
gued for in Section |IV A| and to investigate further the con- 
cern, expressed in Ref. ifTsll . that the ip2 long-range order 
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FIG. 3: (Color online) Order parameter q^^ and heat capacity Cv as 
a function of temperature for J = 1 and V = Q. Left and right panels 
display results for q^^ (L — 2, 3, 4, 5) and for Cv (L = 3, 4, 5), 
respectively. 
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FIG. 4: (Color online) Two plots showing Monte Carlo data for 
J — 1 and D = illustrating the first order nature of the transi- 
tion. The left hand panel is a histogram of the measured energies 
for i = 4 at T/J ~ 0.127, close to the transition temperature. 
The double-peaked structure is evidence for a coexistence region and 
hence an underlying first order transition. The right panel shows the 
peak height of the specific heat, Cv.max, versus the cube of the sys- 
tem size, . The dash line shows a straight line for the hypothetical 
Cv,max oc a + bL'^ in the thermodynamic limit. 



might not survive in the thermodynamic limit, we performed 
Monte Carlo simulations of the nearest neighbour exchange 
only model. Parallel tempering Monte Carlo simulation were 
carried out with J = 1 (and T) ~ Q) for four different system 
sizes, L ^ 2,3,4: and 5, of cubic unit cells of 16 spins. To 
equilibrate the system, 5 x 10^ Monte Carlo sweeps were per- 
formed, followed by the same number of steps to collect data. 
All four system sizes were found to have equilibrated satis- 
factorily according to the criteria discussed in Sec. |llll We 
note that we were unable to obtain well equilibrated results 
for L = 6 even using parallel tempering. 

Figures [3] in and |5] show data from the Monte Carlo simu- 
lations for T>/J — 0. The left hand panel of Fig. [3] shows the 
onset of the ?/;2 order parameter while the right panel shows 
the temperature dependence of the specific heat, Cy, near 
T/ J — 0.127, the estimated transition temperature, for vari- 
ous system sizes. We have found that this estimated transition 
temperature Tc/ J « 0.127 is consistent both with our sublat- 
tice magnetization results (not shown) and those of Ref. ifTsll . 

Both the rate of increase of the heat capacity peak and the 
jump in with increasing L are consistent with a first order 
phase transition in the thermodynamic limit. The left hand 
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FIG. 5: Monte Carlo data for J = 1 and 7? = 0. Histogram of 
the local XY angles for all spins for L = 4 at a temperature T/ J = 
0.1 < Tc. The histogram shows peaks at angles 7i7r/3 illustrating 
that the ip2 states are preferably sampled below Tc/ J ~ 0.127. 



panel of Fig. |4]is a histogram of the measured energies close 
to the transition temperature for L = 4. Its double-peaked 
structure is a clear indication of co-existence and hence of the 
first order nature of the transition. The right panel of Fig. 2] 
shows how the peak height of the specific heat, (7v,max, de- 
pends on the the cube of the system size, L^. For a first order 
transition, one expects Cv.max oc (a + bL^) in the limit of 
large LM The plot illustrates that for L = 4 (L^ = 64) and 
L = 5 (L^ = 125), Cv.max is approaching this expected 
behavior. This provides further evidence for a first order tran- 
sition in this V = (111) pyrochlore XY antiferromagnet. 
Finally, Fig.|5]shows a histogram of the local XY angle aver- 
aged over all spins on all sublattices at T/J = 0.1 for system 
size L = 4. The figure shows six sharp peaks concentrated at 
the tlj2 angles mr/B. We find, in addition, that the spin angle 
on all sublattices are concentrated around one of these angles 
at any given Monte Carlo time. This result therefore demon- 
strates the selection of ip2 states from the continuous mani- 
fold of classical ground states and also that all six magnetic 
domains are sampled in the course of the simulation — a pos- 
sibility facilitated by the use of a parallel tempering algorithm 
in our simulations compared to those of Refs. jlsUISIl . 



D. Further Cases of Order-by-Disorder in Pyrochlores 

We have described in detail the nature of the ground states 
in the (111) XY pyrochlore antiferromagnet and how the clas- 
sical degeneracy is resolved via an order-by-disorder mecha- 
nism. In contrast, the classical Heisenberg model on a py- 
rochlore lattice, which has a much less constrained set of 
ground states, exhibits no phase transition down to zero tem- 
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perature. - To put our results in the broader context of order- 
by-disorder in pyrochlore systems, we summarize in this short 
section two previously studied cases where the degeneracy of 
the Heisenberg antiferromagnet can be lowered by including 
additional interactions, leading to an entropic selection of a 
discrete state of long-range ordered states. 

The first case is the pyrochlore Heisenberg model with both 
isotropic exchange and Dzyaloshinskii-Moriya (DM) interac- 
tionslL32 



= S, • Sj + ^ D„- • (S, X Sj). 



Figure 2 in Ref. Isill gives the D^j which are completely de- 
termined by the lattice symmetry. We have taken a positive 
sign in front of the DM term to denote the so-called indirect 
DM couplings of Ref. Isill . In this case, the single tetrahe- 
dron ground states have four branches. One of these branches 
correponds exactly to the equal (j)a angle ground states of the 
(111) XY model [our Branch 1 of Eq. ©J. The other three 
branches are coplanar spin configurations - the branches are 
distinguished by the three mutually perpendicular normals to 
these planes in the (100) crystallographic directions. The au- 
thors of Ref. Isill observed that Monte Carlo simulations at 
low temperatures lead to the selection of a discrete state of 
states - breaking down the continuous zero temperature clas- 
sical degeneracy down to a Zq symmetry. In the notation 
convention of Ref. lUsll . these are the 1/^3 states. However, 
the nonzero temperature ordering turns out to be more com- 
plicated. As observed very recently in Ref. ||32|l, upon low- 
ering the temperature from the paramagnetic phase, there is 
a phase transition into a V'2 long-range ordered phase fol- 
lowed, at a lower temperature, by the ij^^ ordering reported 
in Ref. jsill . This finding was confirmed within a harmonic 
Holstein-Primakoff computation of the free energy 

The second example of order-by-disorder we discuss is 
the Heisenberg pyrochlore antiferromagnet with both nearest- 
neighbor isotropic exchange and second neighbor interac- 
tions. When the second neighbor exchange is ferromag- 
netic, the finite temperature phase diagram explored by Monte 
Carlo^ exhibits an intermediate phase that appears at interme- 
diate temperatures between the collective paramagnet and the 
low temperature incommensurate multiple-q ordered phase. 
This intermediate phase is partially ordered in the following 
sense. The magnetic structure is layered — each layer exhibit- 
ing collinear spins in one of three perpendicular axes n which 
is common to all the layers. However, the orientation of the 
spins along the n axis is apparently not correlated between the 
layers. The partially ordered phase is selected entropically as 
confirmed by a computation of the free energyi^ The latter 
analytical calculation needs to, and does, include anharmonic 
terms to lowest order because, remarkably, the entropically 
selected intermediate phase is not a local minimum of the har- 
monic approximation to the free energy. 



V. MATERIALS CONTEXT 
A. EraTiaOy 

One motivation for our interest in the local easy plane py- 
rochlore antiferromagnet is the material Er2Ti207. The Curie- 
Weiss temperature 6'cw of this material is negative^ vary- 
ing between — 13 K and — 22 K depending on the temperature 
range of measurement. Er2Ti207 exhibits a transition at about 
Tc = 1.2 Kjii The^ magnetic structure of the ordered phase 
proposed in Ref. \l% has been confirmed by an analysis of 
polarized neutron scattering data.'- 

The spin wave spectrum has been measured recently^ re- 
vealing the presence of an almost gapless mode at q = 
with a linear dispersion for small |q|. This result is consis- 
tent with the magnetic specific heat trend with temperature 
T, '^'''^ The existence of a true Goldstone mode in this material 
is inconsistent with the selection of a discrete set of ordered 
states. The spin wave and specific heat experimental results 
must therefore be providing an upper bound on the size of the 
spin wave excitation gap in this material Specifically, the in- 
elastic neutron scattering spectrum in Ref. 113511 can resolve the 
existence of a gap that is larger than about 1 K and the heat 
capacity measurements of Ref. |39;] show that the heat capac- 
ity varies as down to the base temperature of 450 mK. One 
especially interesting feature of this material is that it appar- 
ently exhibits a field-driven quantum phase transition. Ap- 
plication of a magnetic field in the [110] direction below Tc 
causes a canting of the spins and leads to a zero mode in the 
vicinity of 1.5 T, which has been put forward as evidence for 
a finite field quantum phase transition in this materiali^^i^ 

A greater understanding of the quantum critical point is 
likely to rely on a deeper understanding of the microscopic 
mechanism responsible for the zero field transition in this ma- 
terial. One difficulty with the local XY pyrochlore antifer- 
romagnet as an effective model for Er2Ti207 is that, unlike 
Er2Ti207, the model exhibits a strong first order transition as 
we demonstrated in Section |IVC]'^ Another crucial difficulty 
is that the dipolar coupling of Er'^+ ions is a significant per- 
turbation to the estimated exchange. - The dipolar interaction 
renders the ■4>2 structure energetically unstable and favors a 
distinct ground state, the ?/'4 (Palmer-Chalker) state, as pointed 
out originally in Ref. ifTill . In the next subsection, we concen- 
trate on the dipolar interaction as a perturbation to isotropic 
exchange interactions. Specifically, we consider the possibil- 
ity that thermal order-by-disorder might persist in the pres- 
ence of dipolar perturbations. After that, in Section [TV CI we 
summarize the theoretical constraints that have been placed on 
the problem of understanding Er2Ti207 in zero applied field. 



B. Effect of competing dipolar interaction 

In this section, we look at the effect of introducing pseudo- 
dipolar interactions on the order-by-disorder in the (111) XY 
pyrochlore antiferromagnet. We find that order-by-disorder 
into '4>2 states survives for sufficiently small pseudo-dipolar 
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FIG. 6: (Color online) Figures showing the order parameters q-^,^ 
and q^^ and heat capacities for various values of T)/ J . The top two 
rows, from top to bottom for V/J = 0.5 x lO"'* and V / J = 10"'', 
show q^^ and heat capacity data for i = 3, 4, 5. The third row down 
shows Cv and both g^,, (solid symbols) and q^^ (open symbols) for 
T>/J = 2x 10"'* witli L = 3 (down triangles) and i = 4 (circles) 
data. The bottom row of figures shows q^^ and heat capacity data for 
L = 3 and V/J = 3, 4, 5, 6 x 10"'', with V increasing as indicated 
by the arrows. 



couplings and estimate the maximum value of the dipolar cou- 
pling strength T) that permits 7/^2 ordering. 

The dipolar interaction in the following is taken to act solely 
between nearest neighbors in order to reduce the difficulty 
in equilibrating the system and hence decrease the computer 
simulation time. We study this simplified model because (i) 
in combination with antiferromagnetic exchange, dipolar in- 
teractions between nearest neighbors select the same V'4 mag- 
netic order as the long-range dipolar interactions^" and be- 



cause (ii) we have found, through preliminary simulations of 
Monte Carlo simulations of a model with long-range dipolar 
interaction (not reported), that the general conclusions of this 
subsection below do not depend sensitively on the range of 
these pseudo-dipolar interactions. 

In order to examine the effect of dipolar perturbations on 
the thermal order-by-disorder mechanism, we carried out par- 
allel tempering Monte Carlo simulations of the model with 
Hamiltonian ([T]) with nearest-neighbor dipolar interactions. 
With the introduction of these pseudo-dipolar interactions, the 
simulations were found not to equilibrate sufficiently for sys- 
tem sizes L > 5 for V/ J < 10"'* and poor equilibration was 
found for L > 4 for V/J > 10"''. The results presented in 
this section were obtained from simulations with lO'^ Monte 
Carlo sweeps of which the last 1/lOth was used to compute 
the thermal averages after equilibration. Increasing the length 
of the simulations to 10^ Monte Carlo sweeps was found nei- 
ther to improve the results given below nor to achieve equi- 
libration for larger system sizes than those presented below. 
The simulations were performed on SHARCNET using the 
Saw cluster with Intel Xeon 2.83GHz processors, with the 
L = 5 simulations requiring 10^ Monte Carlo sweeps con- 
suming ~ 17 CPU hours. Each simulation was run on 64 
processors — one processor for each temperature. 

Figure|6]shows the heat capacities and order parameters 
and for different values of V/J (see Section|lIl]for a def- 
inition of g^,2 and (7^4). The two top panels (left and right) 
show the results for V/J = 0.5 x 10""' for L = 3, 4 and 5. 
There is a clear onset of the V'2 order parameter upon lowering 
the temperature and this feature becomes sharper when the 
system size increases. Similar behavior is shown in the sec- 
ond row down in Fig. |6] for a larger value, V/J ~ 10"*. The 
i/'4 order parameter was also measured for these two values of 
V/J and found to be close to zero, within the present finite 
size effects, and displaying no perceptible features around the 
onset temperature of . This is, therefore, compelling evi- 
dence for the thermal order-by-disorder of an ordered ^'2 state 
at nonzero temperature persisting in the presence of pseudo- 
dipolar interactions that break the ground state degeneracy 
in such a way that the thermally selected states are not the 
true ground states. The figures in the third row down (for 
V/J = 2 X 10"'*) reveal strong competition between V'2 and 
tl'4 ordering as indicated by the fact that both order parameters 
are finite but small and suppressed upon increasing the system 
size. We were unable to obtain a clear signature of any long- 
range order in this region of pseudo-dipolar coupling in spite 
of the presence of a heat capacity peak that sharpens with in- 
creasing system size (right-hand panel, third row down). The 
bottom left-hand panel shows the V'4 order parameter for in- 
creasing V/J showing the growing robustness of the transi- 
tion to V'4 order as the pseudo-dipolar coupling is increased. 
For these values of V/J, despite long runs of 10^ Monte Carlo 
steps, the largest system size we were able to equilibrate was 
L = 3. 

The dipolar coupling of nearest neighbor Er"^+ moments in 
EraTisOr is V/{2V2)^ = Am = MoIsjA^b) 74^^^ = 
0.02 K. Whether other (anisotropic) exchange interactions^^ 
acting between the Er^^^ ions are small compared to the 
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isotropic exchange is not known at this time (see Ref. ll40ll 
for a discussion of this problem for the Yb2Ti207 XY py- 
rochlore). 

Nevertheless, a fit to the local susceptibility of Er2Ti207 
in a [110] magnetic field of 1 T (see Ref. 14111 ) has yielded 
an order of magnitude estimate for the exchange coupling 
of about 10"^ K giving (Ain/-'^)cxp ~ 0.2. On the other 
hand, Monte Carlo simulations gives an approximate bound- 
ary between the tp2 and ip4 long-range ordered phases of about 
T)/ J ^ 2 X lO^"' (see Fig.|6l) or, in other words, a condition for 
the appearance of the 7/12 phase of Dnn/ J ^ 0.005. Therefore, 
one is faced with the following conundrum. For the experi- 
mental strength of the dipolar coupling and the estimated ex- 
change coupling, the ordering transition in Er2Ti207 should, 
according to our Monte Carlo simulations, take place into 7/;4 
(Pakner-Chalker) states — but an ordering into 7/12 states is 
observed experimentally. This highlights a problem with the 
thermal order-by-disorder scenario for Er2Ti207. 



C. Towards a model for Er2 Th O7 

It is perhaps worth summarizing the physics problem as- 
sociated with the zero field ordered state of Er2Ti207. The 
long-range ordered magnetic structure that is experimentally 
observed in Er2Ti207 arises via thermal order-by-disorder in 
the local (111) XY antiferromagnet, as discussed in Section 
HV] A large S expansion indicates, furthermore, that there is 
quantum order-by-disorder with the result that the ij}2 states 
are ground states of this model. However, as we have shown 
in Section rVBl the addition of a dipolar interaction (which is 
known to be present and relatively large in Er2Ti207) is in- 
consistent with the selection of ?/'2 ordered states in the clas- 
sical model. In order to resolve this paradox we might be led 
to consider the following possibilities: (i) quantum order-by- 
disorder overcomes the dipolar interaction as conjectured by 
Ref. lITill . (ii) the nearest-neighbor bilinear exchange interac- 
tion is anisotropioi^i^S^ and uniquely energetically selects a 
?/'2 ordered state, (iii) there exist multipolar interactions be- 
tween the Er'^"'" ions accounted for neither in the present work 
nor in previous oneSf^ (iv) the crystal field is responsible for 
the effective anisotropy stabilizing the '4>2 long-range ordered 
statoi^ and (v) that further neighbor interactions might be im- 
portant in this problem. 

It has been shown that allowing for anisotropic exchange in- 
teractions does not lessen the conceptual difficulties in under- 
standing the i}} -2 or dering in Er2Ti207J^ Specifically, as dis- 
cussed in Ref. 111911 . no bilinear nearest-neighbor interactions 
in the idealized (111) XY model on a pyrochlore can bring 
about unique classical zero temperature 7/'2 ground states. The 
11)2 Structure can be obtained, however, via local mean field 
theory by considering the full crystal field spectrum in addi- 
tion to the anisotropic exchange. In this case, the effective 
single ion sixfold anisotropy arises from interaction-induced 
admixing of excited crystal field wavefunctions into the sin- 
gle ion ground state doublet.'- Unfortunately, even in this 
case, the strong long range dipolar interaction between Er''^ 
moments in Er2Ti207 is expected to drive the system into a 



different magnetically ordered (Palmer-Chalker) state, ip^^ 
This is because the energetic selection produced by the dipolar 
interaction inherent to the Er'^+ magnetic moments outweighs 
the six-fold single ion anisotropy effect above (that perturbs 
the XY-like single ion ground state doublet) by at least two 
orders of magnitude. The weakness of the 4'2 selection pro- 
ceeding via the involvement of the excited crystal field states 
is due to the large size of the crystal field gap 10^ K— ) 
compared to the perturbing interactions In short, the com- 
bination of anisotropic exchange and crystal field effects do 
not appear to be able to win against the strong dipolar interac- 
tions and energetically stabilize the experimentally observed 
state rather than the ijj^ (Palmer-Chalker) — state. This 
would seem to rule out possibilities (ii) and (iv). 

In addition, the existence of a ground state crystal field dou- 
blet separated from the first excited states by a gap much 
larger than the scale of the interactions implies that mul- 
tipolar interactions between pairs of angular momenta in 
the microscopic model map to anisotropic bilinear exchange 
couplingsi^ in the low energy effective spin- 1/2 theory <22 thus 
seemingly ruling out multipolar interactions as the micro- 
scopic mechanism behind the zero field long-range order in 
Er2Ti207.^'* Further neighbor interactions (which can also be 
anisotropic) have not yet been investigated. 

We might consider instead a compromise between 
anisotropic exchange and thermal order-by-disorder — that 
important nearest-neighbor exchange anisotropy might more 
than cancel off the nearest-neighbor part of the interaction 
coming from the long-range dipolar interaction. The local XY 
antiferromagnet in the presence of pseudo-dipolar interactions 
with a coupling of the opposite sign [i.e. 2? < in Eq. ([TJ] to 
the physical long-range dipolar interaction preserves a contin- 
uous degeneracy that includes the discrete ?/)2 states;'- How- 
ever, the degeneracy is reduced from the case with only ex- 
change interactions. We might naively expect that the 7/12 
states will generically cease to be the configurations with the 
highest density of soft modes when such anisotropic interac- 
tions are introduced and hence that thermal order-by-disorder 
may not occur, at least not along the lines of the specific mech- 
anism discussed in Section ITV Al However, this naive reason- 
ing may not be correct in general, as we now briefly mention. 

An exception that was noted recently^ is the pyrochlore 
Heisenberg antiferromagnet with Dzyaloshinskii-Moriya in- 
teractions. As mentioned in Section HV Dl Monte Carlo sim- 
ulations and calculations of the free energy for that model re- 
veal that ■02 states are entropically selected at intermediate 
temperatures between the high temperature paramagnet and a 
low temperature coplanar state. The effect of dipolar inter- 
actions was not investigated in Ref. (3^. However, one may 
speculate that in this, or some other model with anisotropic 
exchange, in the presence of dipolar interactions, 7/12 states are 
selected as an intermediate temperature phase which remains 
as a metastable phase upon cooling. In other words, only one 
transition might be observed although two would be observed 
if the system were able to equilibrate. 

The remaining possibility - that the ■4)2 states are the ground 
states of the quantum model despite strong dipolar interac- 
tions classically stabilizing a Palmer-Chalker ijjj^ state - is one 
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that has not yet been investigated. However, in the presence 
of the resulting strong and necessarily nonlinear (e.g. beyond 
1/5) quantum fluctuations, required to overcome the dipole- 
induced ?/'4 ground state, one might expect a substantial re- 
duction of the ordered moment compared to the moment of 
the noninteracting doublet. This expectation contrasts with 
the experimental situation where the ordered moment is about 
3fiB^ compared to a noninteracting moment of 3.8 fiB - a 
mere 20% reduction. Still, such an observation perhaps does 
not rule out quantum order-by-disorder which may therefore 
remains a possible solution to the problem of zero field order- 
ing in Er2Ti207 as originally conjectured in Ref. lITill . We 
hope that a future study will investigate the problem of the 
quantum ground state of the local (111) XY antiferromag- 
net in the presence of dipolar interactions perhaps by carrying 
out an anharmonic 1/5* calculation of the ground state energy. 
Meanwhile, the microscopic mechanism that gives rise to a 
ip2 ordered phase at 1.2 K and with seemingly small quan- 
tum zero point fluctuations in Er2Ti207 remains an open and 
interesting problem in the field of high frustrated magnetism. 

In summary, no mechanism has yet been identified that can 
simultaneously explain how the (dipole-driven) energetic se- 
lection of the ip4 state is avoided in Er2Ti207while allow- 
ing for selection of the iIj2 state below a critical tempera- 
ture Tc ^ 1.2 K. This may indicate that quantum fluctua- 
tions do play a crucial role in the ip2 ordered state observed in 
Er2Ti207 and that a renewed investigation of their effects is 
warranted. 



VI. SUMMARY 

We have presented in Sections and |IV] in some detail, 
the classical ground states, thermal and quantum behavior of 
the pyrochlore (111) XY antiferromagnet with exchange only. 
We have shown that the classical ground states on a single 
tetrahedron have four branches of ground states each with one 
continuous degenerate set of states involving the smooth rota- 
tion of all four spins simultaneously, confirming the previous 
resultji^ From a calculation of the ground states on a single 
tetrahedron, we have inferred the ground states on the py- 
rochlore lattice which include line defects implying that the 
number of ground states scales as L^, where L is the edge 
length of the crystal. 

Monte Carlo simulations of this model confirm that there is 
thermal selection of a discrete set of spin states, denoted ip2 
states,.i2, with ordering wavevector q = from the manifold 
of classical ground states. We have shown, furthermore, that 
this selection occurs to harmonic order in small angular fluc- 
tuations about the classical ground states. In this model, the 
linear spin wave spectrum shows strong similarities with the 
spectrum of eigenvalues of the Hessian. Specifically, the spin 
wave zero modes and the Hessian zero modes appear within 
the same planes in reciprocal space for the '02 states. It fol- 
lows that the quantum zero point energy is minimized at the 
same t/j2 spin configurations that are selected through a ther- 
mal order-by-disorder mechanism, hence confirming the con- 
jecture of Refs. ifiifisli . 



We have considered the effect of including (nearest- 
neighbor) pseudo-dipolar interactions together with the anti- 
ferromagnetic exchange in the classical model at finite tem- 
perature. The (energetically selected) ground states of this 
model are the ?/'4 states'- (also referred to as the Palmer- 
Chalker state>22), given in Section [III so the introduction of 
dipolar interactions produces a competition between energetic 
selection and thermal selection. We have found evidence for 
the persistence of an order-by-disorder transition to a tp2 state 
even when V ^ 0. This finding implies that, in principle, 
a second transition should occur at lower temperatures into 
the '04 (Palmer-Chalker) long-range ordered state since it is 
the classical ground state. However, using Monte Carlo sim- 
ulations, we have found no evidence for such a transition, at 
the very least, because the difficulties of equilibration within 
the ordered ^2 phase prevent the exploration of the space 
of configurations computationally. A similar situation arises 
in a model that tunes between the Heisenberg pyrochlore 
antiferromagnet and the fee Heisenberg antiferromagnet in 
which thermal order-by-disorder, studied using Monte Carlo 
simulations, prevails over the energetically driven ordering 
identified within mean field theory42i No transition was re- 
ported in Ref. ll43ll from the entropically selected ordered 
phase to the energetically favoured magnetic ordered phase. 
In contrast, the high and low temperature phase boundaries 
of an entropically stabilized intermediate (finite temperature) 
state have been identified in the J1/J2 pyrochlore Heisenberg 
antiferromagnet M 

To conclude, we have shown that ■i/j2 long-range order is 
present at low temperatures in the (111) XY pyrochlore an- 
tiferromagnet induced both by thermal and quantum fluctua- 
tions. How such ip2 states, either as a zero temperature ground 
state or as an ordered state at nonzero temperature, arise in 
the presence of (long-range) dipolar interactions remains an 
intriguing question that will require further theoretical inves- 
tigation. 
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Appendix A: Calculation of the Ground States of the XY 
Antiferromagnet 

The zero moment conditions on a single tetrahedron are 

cos (01 ) + COS (02) = COS (03) + COS (04) (Al) 
cos (0'^) + cos (03) = COS (02) + COS (04) (A2) 

cos (0'/) + COS (04) = cos (02 ) + cos (03) , (A3) 
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where (f)'^ = (j)a + and (j)'^ = + ^- The angles 
(j)a are angles in the local coordinate system on sublattice a 
(a — 1,2,3, 4). As discussed in Section |ll] and as previously 
reported in Ref. lfl6ll . the ground states of the (111) XY py- 
rochlore antiferromagnet are characterized by the following 
four set of solutions ("branches"): 



Branch 1 
Branch 2 

Branch 3 
Branch 4 



91 = 03 , 



3 
47r 



(A4) 



In this section, we present a derivation of this result that 
differs from the one in Ref. lfl6ll . Defining Oab = (0a + 0b)/2 
and 5ab = (0a — 06)/2, we first proceed to rewrite the zero 
moment conditions as 

C0s(fTi2) COs((5i2) — COs(cr34) COs((534) = (A5) 

- sin(cri2) sin((5i2) + V3cos(f734) cos((534) = (A6) 
%/3cos((Ti2) sin(5i2) - sin((T34) sin((534) = (A7) 

by using half-angle formulae and then combining the second 
and third zero moment conditions in Eq. ( IA3b to finally obtain 
Eqs. (IA6b and ( IA7b . Our strategy is to eliminate the sums of 
pairs, cr^jy, keeping only the differences of pairs, (5^,y. Thus, 
from Eq. ( IA5b . we get 



sin^(cri2) = 1 



COs(cr34)cOs((534) 



COs((5i2) 



Substitute into Eq. ( IA6b to get 



(1 - sin2((5i2) - cos2((534) cos2(cr34)) sin2(5i2) 

= 3cos2(a34)sin2(534) (l - sin2(,5i2)) (A8) 

Then, squaring Eq. (IA6b and Eq. (IA7b . and adding the result, 
we obtain 

sin2(5i2) = Q sin2(a34) + Zco^^{g^S\ sin'(<534) (A9) 



which we can otherwise write as 



cos (0-34) = - 



3 /sin2(,5i2) 1 



sin2(^34) 



(AlO) 



So now we can proceed with what we set out to do; substi- 
tute Eq. ( lAlOb into Eq. (IA8b leaving us, after re-arranging and 
cancelling off a 003^(^12) term, with 

sin^((5i2)sin^((534) = 
1 



[3sin2(5i2) - sin2(534)] [3sin2(<534) + 



(All) 



It follows from this last equation that 

^(sin2(^34)-sin2(^i2)) =0. 



Hence, the most general form for the set of 0a angles is 

(01,02,03,04) = (0 + 6',0,V'±6',V'). 

Substituting this into our original zero moment formula, 
Eq. ( IA5b . which we write here again: 

cos (0i) + cos (02) = COS (03) + cos (04) , 

shows that one must have either 

(01,02,03,04) = (0 + 6',0, + 6',0) 



or 



(01, 02, 03, 04) = (0 + 61, 0, 0, + d) 

if 6* is nonvanishing or 

(01,02,03,04) = (0, 0, ±0, ±0) 

if 6* = 0. Thus, the angles must occur in pairs. One can 
now return to original zero moment conditions, Eq. ( IA3I ). and, 
with the knowledge that the angles must occur in pairs, obtain 
the four branches of ground states in Eq. ( IA4I ). For example, 
suppose that the pairs of angles occur in the configuration 

(?/;,0, 0, ^) 

then Eq. (IA3I ) gives the branch 



47r 47r 
T ~ ~3~ ~ 



Appendix B: Brillouin Zone Integration 

In the calculation of the free energy in Eq. ([8]l, in particu- 
lar the last term of that equation, we take the Brillouin zone 
sum over to an integral as applicable to the case of an infinite 
lattice. In general, we expect 

where the integral is taken over the Brillouin zone of volume 
i^BZ and A^p is the number of primitive cells. We thus have 



E 



4(4(2^)3) JBZ 



The edge length of the cubic unit cell, a, has been set equal 
to one. The two factors of one quarter come about because 
(i) 4Ap = N where N is the number of spins and (ii) the 
Brillouin zone volume is 4(27r/a)3. 
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Appendix C: Equipartition argument 

In this section, we revisit the Monte Carlo simulations pre- 
sented in Ref. fisll . In particular, we consider Fig.|6]in that 
work showing the average energy in the ordered phase at low 
temperatures as a function of the system size. The authors of 
Ref. in found that 

E 1 

We find that the existence of chains of zero modes in the (111) 
XY pyrochlore antiferromagnet is sufficient to constrain the 
coefficients a and /?. 

Let us consider a cubic cell with edge length L with each 
cubic unit cell of unit edge length. The number of spins is 
N = IQL^. Consider a square face of a single cubic cell. 
There are two chains beginning on the face that alternate be- 
tween sublattices a and h. The number of such chains passing 
through the sample is iVchains = 2L^. There is one degree of 
freedom per spin so the average energy would be {l/2)NkT 



were it not for the fact that the spectrum of modes about the 
ip2 States has planes of zero modes. We shall assume that the 
zero modes (in the harmonic spectrum) are actually resolved 
by a quartic contribution to the energy when looking at the 
higher order corrections to the Hamiltonian. The zero modes 
correspond to rotations along two classes of ab chains so the 
number of such modes is 2A^chains. Thus, the average energy 
is 

-r^ = liN- 2A^chai„s) + 7 (2iVchains) = Ja^^ ^A^chains- 

KbJ 2 4 2 2 

Then, because A'^chains = SL^ = 2(iV/16)2/3, we obtain 

NkBT 2 162/3 2 16 L 

The coefficient in front of 1/L from Monte Carlo simulations^ 
is approximately 0.0636 which is in good agreement with the 
calculated 1/16 = 0.06250 value above. 
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In more detail, suppose that the microscopic interactions between 
magnetic ions are two-body interactions that respect the lattice 
symmetries. Suppose, furthermore, that the crystal field spectrum 
has a ground state doublet and a gap to excited single ion states, 
A, that is much larger than the energy scale associated with the 
interactions. Finally, we suppose that we are interested in temper- 
atures T <C A. Then the excited crystal field can be neglected and 
the effective theory at low energies is the projection of the micro- 
scopic Hamiltonian onto the product of crystal field doublets on 
each magnetic ion. Since the space of states at each magnetic site 
is two dimensional, all site operators are linear combinations of 
Pauli spin operators regardless of the nature of the bare micro- 
scopic interactions. Therefore multipolar interactions project onto 
bilinear exchange couplings. 



